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Abstract 


A  perturbation  approach  to  the  problem  of  determining 
the  electromagnetic  field  reflected  from  the  ionosphere  is 
presented.  It  is  shown  that  by  treating  the  earth's  field 
Of  magnetic  induction  as  a  perturbation^  Maxwell's  equations 
and  the  constitutive  relationship  are  considerably  simplified. 

First  order  approximations  for  the  field  reflected  from 
a  homogeneous  sharply  bounded  ionosphere  are  found. 

The  perturbation  method  to  be  developed  has  certain 
limitations  which  are  discussed. 


I.  INTRODUCTION 

I, I.  The  Reflection  Problem  formulated 

The  problem  of  electromagnetic  reflection  from  the 
ionosphere  will  now  be  formulated  in  very  general  terms* 

The  coordinate  system  is  arranged  such  that  the  z-axis 
is  vertical  and  the  earth's  field  of  magnetic  induction 
lies  in  the  x-z  plane.  The  ionosphere  is  assumed  to  be 
horizontally  stratified,  non-magnetic ,  and  confined  to  the 
half  space  z>0.  The  region  z<0  is  assumed  to  be  free 
space.  It  is  also  assumed  that  a  known  electromagnetic 
source  is  located  in  the  free  space  region  z<0.  let  the 
source,  an-d  therefore  the  field  quantities,  have  a  harmonic 
time  variation  of  the  form  Since  the  ionosphere  has 

been  assumed  horizontally  stratified,  the  electron  density 
and  collision  frequency  will,  in  general,  be  a  function  of 
the  coordinate  z.  These  functions  are  assumed  to  be  at  leaet 
piece-wise  continuous. 

For  the  purposes  of  ionospheric  wave  propagation.  Max¬ 
well's  Equations  (MKS  units)  and  the  constitutive  relation^ 
ship  may  be  written^ 


V  X  B  ■  iooji^H 

V  X  H  “  i«6  (5  ♦  T  )  *  7 
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where 

N  =  Electron  density 
e  «  Electronic  charge 
in  =  Electronic  mass 

"I  =  Earth's  field  of  magnetic  induction 
V  =  Collision  frequency 
j  =  Current  density  of  the  source 
a  »  Angular  frequency  . 

The  plane  z  »  0  is  an  interface  between  two  media;  namely 
the  ionosphere  and  free  space.  If  j  =  0  on  the  plane  z  =  0, 
then  the  tangential  components  of  E  and  H  must  be  continuous 
on  the  plane  z  «  0. 

The  problem  statement  is  extremely  simple.  Given  the 
current  density  of  the  source  and  an  ionospheric  model,  find 
E  and  H  such  that  equations  (l)  and  (2)  are  satisfied  and  the 
tangential  components  of  E  and  H  are  continuous  on  the  plane 
z  -  0, 

I.JI  General  Remarks 

Witltout  loss  of  generality,  the  source  may  be  replaced 

by  a  point  dipole.  The  dipole  field  may  be  represented  as  an 

2  3 

integral  over  plane  waves  *  . 

A  possible  method  of  attack  would  be  to  find  the  field 
reflected  from  the  ionosphere  when  a  plane  wave  is  incident 
at  an  arbitrary  angle  upon  the  ionosphere.  Then,  by  integra¬ 
tion  over  the  incident  angles,  the  reflected  field  of  a  point 


If 


dipole  may  be  obtained.  But  first  a  plane  wave  solution  must 
be  found. 

I. Ill  Problem  Obtaining  a  Plane  Wave  Solution 

A  great  deal  of  work  has  been  done  with  plane  waves  in 
the  ionosphere.  However^  most  methods  are  not  of  sufficient 
generality  to  be  of  use  in  the  solution  of  this  problem.  It 
seems  three  basic  methods  of  finding  solutions  for  a  plane 
wave  of  arbitrary  incidence  are  available;  other  methods  are 
seemingly  special  cases  of  these  three.  These  methods  are 
listed  below: 

1)  Booker  Quartic^ 

The  ionosphere  is  first  assumed  homogeneous.  The  roots 
of  the  Booker  Quartlc  give  the  propagation  constants.  Two  of 
the  roots  are  aesoclated  with  the  up-going  plane  waves  and  the 
other  two  roots  are  associated  with  down-going  plane  waves. 

'The  solution  is  extended  to  an  inhomogeneous  ionosphere  IB^ 
using  the  W.K.B.  approximating  technique. 

2 )  First  order  coupled  equations 

Clemmow  and  Heading  present  a  set  of  first  order  coupled 

equations  for  the  quantities  E  ^  E  ^  H  and  H  •  The  equations 

X  y  X  y 

t 

uncouple  in  a  homogeneous  medium.  Keller-  gives  a  detailed 
method  of  solving  these  equations. 

3)  Integral  equation  variational  technique^ 

A  variational  technique  is  presented  for  the  reflection 
and  transmission  matrices.  The  power  of  this  technique  is 


-u- 


that  the  detailed  solution  of  the  fields  In  the  ionosphere 
never  need  be  foundb 


The  methods  mentioned  above  have  the  disadvantage  of 
being  rather  laborious.  Also  the  solutions  obtained  by- 
these  methods  involved  the  angles  of  incident  in  a  compli¬ 
cated  fashion.  Hence,  the  integrals  over  the  resulting  plane 
waves  cannot  be  evaluated  short  of  numerical  means. 


* 


II.  DEVELOPMENT  OF  PERTURBATION  METHOD 


Maxwell's  Equations  and  the  constitutive  relationship 
are  written 


Vx  E 


icoji  H 
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Vx  H  +  icoe  (B  ♦  I  )  *  j 
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Treat  the  parameter  Y  as  a  perturbation  parameter. 


Since  Y 


is  directly  proportional  to  the  earth's  field  of 


magnetic  induction 


equivalent  to  treating  B 


B  ,  treating  Y  as  a  perturbation  is 


as  a  perturbation. 


Expanded  E,  H,  and  P  is  a  power  series  in 


YJ  about 
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The  expansions  (H)  are  substituted  into  equations  (l), 
(2),  and  (3)  and  the  vector  coefficients  of  like  powers  of 
I?  I  are  equated  to  obtain  the  following  relations: 
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Vx  H"  +  icoe  +  -7  )  =  6 

0  6  no*^ 
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n  >  0 


U  7^ 
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22-  -:^|--U^e'‘-  o(of‘E”"^)-  lUcr  X  e“"^  n>  2 

^  ’■o  ^  ‘o  (9) 


where  .  « 

^  1^1  _  1  n  =  0 

T  ^  Y  cr  6  = 

'  0  n  /  0  . 

Equations  v?),  (8),  and  (9)  are  used  to  eliminate  from 
equation  (6).  Equation  (6)  becomes,  with  nl^Os 

V  X  *  icoe  E^  -  l" 
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n  >  2  . 


The  problem  of  finding  the  solutions  of  equations  (l) 
and  (2)  has  been  reduced  to  finding  the  solution  of  the  follow 
ing  infinite  set  of  equations: 


“7 


Vx  b“  =  (lO) 

Vx  1*^  +  =  T“  (11) 

The  following  observations  are  made: 

1)  If  is  known,  then  equation  (lO)  may  be  used  to 

determine  Hence,  only  need  be  found* 

2)  and  are  bounded  since  ^  and  H  are  bounded.  If 

j  ■  1°  =  0  on  the  plane  z  =  0,  then  the  tangential  com¬ 
ponents  of  and  (n>0)  must  be  continuous  on  the 

plane  z  =  0.  The  validity  of  this  statement  may  be 
deduced  from  equations  (lO)  and  (ll)  in  the  usual  manner* 

3)  l”  depends  on  and  E*^”^*  Therefore,  if  and 

are  known,  is  known* 

U)  The  equations  (lO)  and  (ll),  which  determine  take  the 

form  of  Maxwell’s  equation  in  an  isotropic  medium  with  a 
source  current  density 

Equations  (lO)  and  (ll)  with  n  =  0  are  used  to  find  ¥°, 
Once  E°  is  found,  I^  is  known*  Then  equations  (lO)  and  (ll) 
are  again  used  but  with  n  ■  1  to  find  E^.  This  process  can. 
be  continued  to  calculate  as  many  E^'s  as  desired. 

In  summary,  the  problem  of  finding  the  solutions  of 
Maxwell’s  equations  in  an  anisotropic  medium  is  reduced  to 
finding  the  solution  ef  an  infinite  set  of  equations  by 
treating  |Y  as  a  perturbation*  The  n-th  equation  of  the  set 
takes  the  lorm  of  Maxwell’s  equations  in  an  isotropic  medium. 
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When  and  have  been  determined,  becomes  known, 

and  E^  can  be  found,  is  subject  to  the  same  boundary 

conditions  as  E. 

The  physical  interpretation  of  the  perturbation  series 


”e^  is  interesting.  E°  is  the  solution  of  the  pro- 


n^O 

pagation  problem  with  no  earth’s  magnetic  field.  The  terms 
“1  —  2  2 

YE,  Y  E  ,  etc,  are  correction  terms. 


II. I  Limitations  of  the  Perturbation  Approach 

The  convergence  properties  of  the  series  (ii)  cannot  be 
known  until  the  general  term  is  found.  In  most  practical 
applications  the  determination  of  the  general  term  would  be 
at  best  laborious.  However,  if  j  y|«1  the  first  few  terms  in 
the  series  should  be  sufficient  to  represent  the  solution 
accurately.  The  charge  to  mass  ratio  of  the  electron  and  the 
earth's  field  of  magnetic  induction  are  fixed  quantities. 
Therefore,  the  only  parameter  in  j ij  that  may  be  varied  is  w. 
Thus,  j  Y  1  implies  high  frequencies. 

As  has  been  previously  mentioned,  equations  (lO)  and  (ll) 


take  the  form  of  Maxwell's  equation  in  an  isotropic  medium. 
However,  since  p.  is  a  function  of  position  the  medium  is  inho¬ 
mogeneous,  Thus,  the  positional  dependance  of  p  must  be  such 
that  Maxwell's  equations  can  be  solved. 


III.  PERTURBATION  THEORY  APPLIED  TO  A  S PECIaL  C A S E 


The  ideas  developed  in  the  last  section  will  now  be  used 
to  find  first  order  correction  terms  for  the  field  reflected 
from  a  sharply  bounded  homogeneous  ionosphere.  Let  the 
interface  between  the  ionosphere  and  free  space  be  the  plane 
z  =  0,  A  vertical  point  dipole  is  located  at  (0,  0,  -d) , 

The  earth's  field  of  magnetic  induction  is  assumed  to  be 
vertical.  Since  the  perturbation  approximation  is  applicable 
only  to  high  frequencies,  collisions  will  be  neglected. 

Equations  (lO)  and  (ll).  With  n  =  0,  are 


where 


V  X  E°  = 

V X  =-ioo6^|i^E°  +  j 
2  ('I  -  X-i  z  >  0 

“  ll  J  ZCO 


(12) 

(13) 


6(r  -  r)  =  Dirac  delta  function. 


1  ®  Unit  vector  in  the  z-direction  . 

*z 

In  general,  the  symbol  Ij^  will  mean  a  unit  vector  in  the 
direction* 

If  11°  is  eliminated  from  equations  (l.)  and  (l  )  the 
result  is 

Vx  (Vx  1°)  -  k^l°  (U) 
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2  2 


let 


T 


z  >  0 
z  <  0 


Then  it  follows  from  equation  {ih)  that  E*  and  E^must  satisfy 


?x  (Vx  E*)  -  k^E^  =  0  (15) 

Vx(Vx^”)  -k^E  =  itoti  j  .  ( l6) 

o  o 

Let  us  introduce  two  new  vectors,  n*  and  n  ,  defined  by 


=V(V*«*)  ♦  ,  (17) 

IJ*  =  V(V‘n')  ♦  k^n"  .  (18) 

Then,  from  equations  (15)  and  (l6). 


(V^  *  k^)  TT* 

.  k^)  n- 

o 


(19) 

(20) 


z 


Since  j  has  only  a  z-component. 


n 


t 


0, 


n*  =  0,  and 

y  ’ 


From  equations  (l7),  (18),  and  (12)  the  following  express¬ 
ions  are  found  for  the  field  components  in  cylindrical  coordin¬ 


ates: 

E*  ^  i_  tJL* 
r  3r  3z 


E 
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3  3n 
3r  Fz 


4!  'f* 

-*^18  8n 
<p  r  d<p  dz 


13^ 
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H  = 
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The  conditions  E 


and  Eg  =  Eg  evaluated  at  z  =  0 


imply 


_9  ^ 
9r  9z 


-i.  ill 

9f  9z 


9  9  n  _  9  9  n 
9  <p  d  z  9  (p  9  z 


z  =  0 


z  =  0  , 


(21) 

(22) 


Equations  (2l)  and  (22)  may  be  integrated  with  respect  to 

r  and  <p;  respectively.  In  both  cases  the  constants  of  integra- 

9  Jv* 

tion  must  be  zero  since  ^  tends  to  zero  as  r  tends  to  infinity. 


In  a  similar  manner  it  is  found  that  the  conditions 


H*, 

r*  9 


H 


r'9 


evaluated  at  z  =  0  are  satisfied  if 


2  ♦ 

n  R 


Thus  it  can  be  concluded  that  if  n"  are  such  that 


i!L 

9z 


dz 


2  * 

■  R 


(23) 

(2h) 


then  the  tangential  components  ^  and  "H  will  be  contibuous  on 
the  plane  z  =  0. 

T  4.  -  P  4.  S 

Let  n  ■  R  ♦  R 


where 


iwii^I 

(V^  ♦  kf)  e( 


r„  -  r) 
o 


(25) 
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(V^  k.)  Tt  =  0  , 


(26) 


n 

The  solution  of  equation  (25)  ‘  is 


P 

'  '  r~ 

Uii 


ik  r  -  r 
0  o 


1-  ‘  --ol 


where  r  =  -  1  d  • 
0  z 


ik. 


r  -  r 
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If  the  Sommerfeld  integral  representation  for 


r 

is  used,  ft  may  be  written: 

z  =  p  r  j  ^ 


r  -  r. 


(27) 


where 


lin 


J^(Xr)  =  Cylindrical  Bessel  function  of  the  order 


Equations  (26)  and  (l9)  are  satisfied  if 

e  °°  Q  "Y I 

^  P[  f^  (X)J  (Xr)e  XdX  z^  0 

Jo  ° 


(28) 


-Yd 


i*  ^  P  r  f*  (X)e  J  (Xr)e  *'  XdX  z>  0 

Jo  ° 


where 


(29) 


s  * 

The  functions  f- (x)  and  f  (x)  must  be  determined  from  the 


boundary  conditions. 

The  following  points  should  be  noted: 
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p 

1)  The  function  tt  is  singular  at  (O,  0,  ^  d)  since  the 

integral  (26)  falls  to  converge  at  this  point.  This 

singularity  is  caused  by  the  presence  of  the  dipole, 
v  S  ’*■ 

2)  The  functions  ti  and  n  are  continuous  throughout  the 
regions  z^O  and  z  >  0*respectively . 

3)  All  three  functions  obey  the  Sommerf eld  radiation  condi¬ 
tion  at  infinity. 


By  applying  the  boundary  conditions  expressed  by 
equations  (23)  and  (2li)  it  is  found  that 


f^(\)  =  - 

H  Y  ♦  Yg 

,3a,  -  i  - 


(30) 


(31) 


These  solutions  were  first  obtained  by  Sommerf eld  in  his 

7 

1909  paper  . 

Let  F(X)  -  Pf^d) 


n 


F(X)J  (Xr)e 

P 


dX 


z  >  0  . 


Equations  (lO)  and  (ll)  with  n  =  1  are 


V  X  E^  = 

(31) 

Vx  ^  -  lose  11^#  ♦ 

• 

(32) 

Eliminating  from  equations  (3l)  and  (32)  gives  the 


result 
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Vx  (Vx  =  i£^  (33) 

Vx  (Vx  1^')  -  k^f^“  -  0  ,  i3h) 

where  . 

,  E^  z  >  0 

®  ■  ^1-  z<  0  . 

We  now  introduce  the  vectors  IT^*  defined  by 


E^“  =  V(V*«^"')  *  ^o"^” 

into  equations  (33)  and  (3li)  to  obtain  the  following  wave 
equations: 


(V^  ♦  k^) 

-1+  _  ^^0  jl 

n  _  -  x  ■  ■  1 

(35) 

(V^  ♦  k2) 

=  0 

(36) 

=  CO  e 

X  0  X  E° 

Cf  =  1 

or  X  E  =  1  X  (l  E*  ♦  1,2*)  =  1_  X  I  E* 
z  rr  zz  z  rr 

^  X  E  =  1  B* 
cp  r 


Thus  equation  (35)  becomes 
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where 
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(38) 

(39) 
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The  integral  in  equation  (39)  may  be  evaluated  by  using 
8  9 

the  following  facts  *  : 

iklr-r^l  oo 

1,  ^ - -  -  y  £  cos  m  (<p-(p^)  r  J  (Xr)J  (Xr^)  x 

—  —1  m  [  tn  m 
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e  _ XdX 

Ye 


where 


m 


ri\  m  =  0 

V  2  j  in  0  . 


2. 


j: 


1  l^_l  -  6(X-X^) 


)J  (X'^r'^)r'‘'  dr 


m 


3.  1  ^  1  sin  <p  ♦  1  cos  cp  . 

^  X  y 


-IP  iX 

1C 


ly  r  F(x)  J.(Xr)q(z)e  ®  XdX 
2  Jo 


(Ui) 


where 


q(0  •  ^  ♦  2:^ 


,1P  .  1  « 

q>  <P 


2  H 


J  j  F(X)  (Xr)  q  (z)  e 


XdX  .  (Ij2) 
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It  will  be  found  that  in  order  to  satisfy  the  boundary 
“IS  -1- 

conditions,  n  and  ti  Will  have  to  have  a  z-component  as 
Well  as  a  (p-component 


-1*  IP  ^  ^IS  ISx 

"  ■  (0.  %  *  %  »  *2  >  . 


il- 


The  tangential  components  of  the  fields  arei 


-li 

n 

(li3) 

<p 

1  ) 
r  36 

CUU) 

«r*  = 

,  2 

/  ^  an^"  3n^-  \ 

k 

1  . 9  . 

(li$) 

l“Vo 

1  r  3<p  3z  / 

._ki 

(U6) 

9 

icap. 

n  * 

1+  1- 

The  condition  E  *  E  evaluated  at  z  r  0  implies 
r  r 

■  ^v.  2-0  .  (47) 

Equation  (U?)  is  integrated  with  respect  to  r  to  obtain 

V*n^*  «  z  ^  0  .  (Ii8> 

The  constant  of  integration  is  zero  since  tends  tP  zero 

as  r  tends  to  infinity. 

By  using  the  continuity  of  E^,  equations  (UU)  and  (U8), 
it  is  found  that 

’  Z  0  .  (It9) 
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T^e  continuity  of  and  equation  (it6)  implies 


2  1+  1- 
"'z  "  "z  » 


z  =  0 


(50) 


Finally,  by  usi-ng  equations  (U5)  and  (50)  and  the 


continuity  of  ,  it  is  found  that 
i+  ^  1- 

«  an"  an 

.2  _ <p  _  (p 


a  z 


a  z 


z  •  0 


(51) 


The  following  two  boundary  conditions  are  used  to 
determine  the  arbitrary  functions  that  will  appear  in  the 
integral  representation  of 


2  1+ 
(p 


2 


,1- 

<P 


an 


1- 

■<p 


(52) 


az 


az 


When  n^  has  been  found,  the  two  boundary  conditions 


v.-ji*  .  V-?1- 


2rl+  =  ,1- 

II  n  *  n 

~  tf  rw 


(53) 


are  used  to  find  the  arbitrary  functions  that  will  appear  in 

1± 

the  integral  representations  of  n  • 

It  may  be  verified  by  direct  substitution  that  the 
relations 


.  «  f  oScx)J,(Xr) 

2l»^  Jo  ^ 

oo  yZ 

^  j  G“(\)J^(Xr)  e  XdX,  |d|  >  |z|  z<0 
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00 


00 


$"(X)J^(Xr) 


e  XdX,  z>  0 


z  <"  0 


satisfy  equations  (38)  and  (36)  if  G  (x)^  G  (X),  #  (X),  and 
#  (X)  are  integrable  functions. 

The  boundary  conditions  ($2)  and  (53)  are  used  to  find: 


.  + 

#  =  $ 


G‘(X)  - 


=  0 

tt^F(x) 

Y+Ye 


G 


+ 


1 

Y+Ye 


Since  the  fields  in  the  region  |  dj  ^ 


z<0  are  the 


observable  quantities,  the  fields  in  the  region  z>0  will  not 

be  considered  further,  except  to  mention  that  lim  =  •  . 

X^l 

—1+ 

However,  lim  E  is  a  finite  quantity. 

I^l 


ni  -  iXP 


L 


®  J^(Xr)  e  "Y(d-z) 


0  (ji^  Y'»'Ye)^Y*Ye) 


X'^  dX,  (  d  >  z  )  (z<0).  (5U) 


Since  =  1  and  is  not  a  function  of  q)j 

<p  <p  <p 

— 1—  2  1  — 

B  ^  k  1  ,  Thus  a  first  order  correction  term  for  the 

o  <p  (p 

reflected  field  has  been  found. 

The  integral  (51i)  must  now  be  evaluated. 

Let  d  -  z  ^  , 
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Q(\)  *  — g  ■  - 

Ci»  v+YgHy+Yg) 

Since  Q(X)  is  an  even  function  of  X  and 

2J^(Xr)  *  (Xr)  +  (-Xr),  integral  (5U)  may  be 

written 

Q(X)  (Xr)  e  dX  (55) 

where  (Xr)  is  the  Kankel  function  of  the  first  kind. 

The  integral  is  not  yet  uniquely  determined  because  there 
are  branch  points  at  X  =  ±k^  and  X  =  ±k.  The  integrand  is 
made  unique  by  introducing  branch  cuts.  These  cuts  and  the 
path  of  integration  are  shown  in  Fig,  I  under  the  assumption 
that  (1  <  1. 

With  the  plane  cut  as  shown,  the  poles  of  the  function 
Q(X)  lie  on  the  second  Riemann  sheet. 

In  general  the  integral  (55)  is  difficult  to  evaluate. 
Therefore,  the  problem  will  be  simplified  further  by  assuming: 

1)  Since  the  reflected  field  is  of  interest  when 

z  is  comparable  to  d,  and  d  is  very  large  as  compared 
to  a  free  space  wave-length  (perturbation  theory  is 
only  applicable  to  high  frequencies),  ^  is  large  and 
positive, 

2)  Since  reflections  occur  in  that  region  of  the 
ionosphere  where  the  electron  density  is  such  that 

X  ^  1,  the  limit  of  equation  (55)  as  X  tends  to  unity 
will  be  taken.  It  is  assumed  that  the  limit  may  be 
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taken  inside  the  integral# 


Then  equation  (55)  becomes 


aj"  =  ^  f  Q^(X)M^^^  (Xr)  e  ^^dX 
^  J  — QO 


(56) 


where 


Q  (X)  «  lira  Q(X) 
X-»l 


X 

X  ♦  Y 


By  using  the  asymptotic  form  of  the  Hahkel  function  and 

the  substitution  X  «  k  sin  9,  the  integral  (56)  is  trans- 

0 


formed  to 


itt  n/2-i® 


,1„6 


ik^[r  sinS+gcosoj 


cos  6  e 


'k^  r  sin  9 
0 


(57) 


This  integral  may  be  approximated  by  using  the  saddle 
point  method.  To  this  end,  let  r  =  R  sin  a  and  ^  ■  R  cos  a 
to  transform  (5?)  to 


^  /T  e'  f  IULL 

2  ^  «  i  VTT 


Q  -i0  ik  R  cos  (0-a) 

^  e  ®  d0 

sin  e  (53) 


where  the  contour  C  begins  at  -  n/2  ♦  a  +  i®  and  ends  at 
n/2  *  a  -  i®  . 


®(0  )  ■  ik  R  cos  (9  -  a) 
o 


<p'(9g)  =  -  ik^  R  sin  (9g  -  a)  -  0 
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Therefore,  there  is  a  saddle  point  at  0  =  0g  =  a. 

0*(e)  =  -  ik  R  cos  (0  -  a)  . 

In  the  neighborhood  of  the  saddle  point 


<p(0  ) 


ik  R 
0 


ik  R 
o 

2 


(59) 


(6  -  a)  =  s  e^'* 
ik  R 

<p(0 )  »  ik^R  ^  °  (cos  2  V  -  i  sin  2  v) 

Im  <p(0 )  =  k  E  -  — s  cos  2  V 


The  contour  C  must  be  deformed  into  a  contour  c'  which 

has  the  property  that  along  C'  Im  <p(0)  *  Im  q>(a).  In  the 

neighborhood  of  the  saddle  point  this  condition  requirss 

V  =  ±  n/k*  The  sign  of  v  is  chosen  such  that  R  <p(0 )  <  0  when 

e 

s  >0,  In  this  way  it  is  found  that  v  •*  -  The  contour 

C'is  shown  in  Fig,  II. 

From  equation  (59)  equation  (58)  becomes 

k  R  ,  . 

,  p  _  in  ,  ik  R  --S 

_1-  ,  ^0“  /T"  Tf  f  sin  6  cos  0  e”^®  e  °  e  e  ^  ds  . 

♦  “  i,  Vie/  3ln  e 


Hence,  to  a  good  approximation^ 


sjn  g  cos  a  e 

Vk  r  sin  a 
»  ■  0 


ik  R  - 
e  ^  e 


k  R 
0 


2 

ds 


in 

P  e  ^  cos 


-io 

9  e 


ik  R 
0 

e _ _ 

R 
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The  use  of  the  asymptotic  form  of  the  Hankel  function 

is  justified  if  k  r  sin  a»l.  Because  of  the  way  in  which 

0 

the  angle  a  was  defined,  it  is  seen  that  takes  the  form 

of  a  Spherical  wave  centered  about  the  dipole  image  located 
at  (O,  0,  ♦  d)  . 

In  summary,  the  expression  (60)  holds  if; 


1.  X  =  1 

2,  R>>1  and  k  r  sin  a»l. 

o 

Since  lim  f^(X)  =  -  — 

X-1  ^ 


Therefore,  in  the  region 


d  < 


2  <  0  ,  the  total 


reflected  n  vector  may  be  written 


-  P 


1  * 
z  (p 


Y  i  cos  a  e 


-ia 


]t 


ik  R 
0 


(61) 


The  reflected  field  intensity  B  may  be  obtained  by 
operating  on  by  V(V*)  +  k^  , 
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IV.  POLARIZATION  OF  THE  H  VECTOR 


The  polarization  of  the  reflected  field  is  of  some 
interest.  In  what  follows  it  proves  convenient  to  introduce 
Standard  spherical  coordinates.  This  change  is  easily  made 
in  equation  (6l)  by  replacing  a  by  n-^0  where  0  is  the  co¬ 
latitude.  With  these  changes  equation  (6l)  becomes 

i(e  +  n/2) 


n  ^  -  pfl  ♦  1„  Y 
L  z  <p 


ik  R 
o 


cos  0j  ^  .  (62) 


It  should  be  noted  that  these  angles  are  measured  with 
respect  to  a  coordinate  system  whose  origin  is  at  the  image 
point  (0,  0,  *  d) . 

Since  1  =  1_  cos  0  -  1*  sin  0,  the  components  of  the 

z  R  0  ' 

reflected  IT  vector  are 


ik  R 
o 


‘R 


-  P  cos  6 


Kg  =  P  sin  0 
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ik  R 
o 


<P 


-  P 


i(e  +  n/2) 


ik  R 
.  o 


cos  0 


Using  the  definition  of  the  quantity  P  =  — — ^ 


I  and 


iov 


y  X  n 


H  the  components  of  H  in  the  radiation  region 


are  found  to  be 


ik  I 

^0  -  nih  '  ®  « 


iC0  *  n/2) 
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(p 


ik  I 

■  sin  0 


ik  R 
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Thus 


-icot 
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k  I  I 

jY|  Gos  6  cos  (k^R  +  0  +  n  -  cot) 


e  (p 


-icot 


k  I 

1-°=-  sin  6  cos  (k  R  +  ^ 
UnR  0  2 


cot) 


Let  be  the  angle  between  -  Ig  and  the  major  axes  of  the 
polarization  ellipse^ 


After  a  simple  calculation  it  is  found  that 
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tan  2\|c 


2  1 
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sin^0 

cos  0 
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cos  0 

-  sin^0 
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V.  SUMMARY 

It  has  been  shown  that  by  treating  the  earth's  magnetic 
field  of  induction  as  a  perturbation.  Maxwell's  equations  and 
the  constitutive  relation  are  simplified.  The  partial 
differential  equations  which  determine  the  n'th  coefficient 
of  take  the  form  of  Maxwell's  equations  in  an  isotropic 

medium.  However,  the  medium  is  still  inhomogeneous.  Since, 
in  general,  it  is  not  practicable  to  calculate  correction 
terms  of  order  greater  than  the  first,  the  frequency  must  be 
high . 

A  first  order  correction  term  has  been  found  for  the  H 
vector  reflected  from  a  sharply  bounded  homogeneous  ionosphere. 
The  earth's  magnetic  field  was  assumed  to  be  vertical  and  the 
source  was  taken  to  be  a  vertical  point  dipole  located  in  free 
space  a  distance  d  below  the  interface. 

The  correction  term  is  expressed  in  terms  of  a  Fourier- 
Bessel  integral.  In  order  to  evaluate  this  integral,  the 
limit  as  X  tends  to  unity  was  taken  and  the  saddle  point 
method  was  applied.  Finally  the  angle  between  -  1.  and  the 
major  axes  of  the  polarization  ellipse  of  the  reflected  H 
vector  was  stated. 
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